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^ ■ Abstract 

We consider a hierarchy of many particle systems on the line with polynomial 
(f-) ' potentials separable in parabolic coordinates. Using the bax representation, writ- 

ten in terms of 2 x 2 matrices for the whole hierarchy, we construct the associated 
linear r-matrix algebra with the r-matrix dependent on the dynamical variables. 
A dynamical Yang-Baxter equation is discussed. 



1 Introduction 



It is known that it is possible to provide a 2 x 2 Lax operator satisfying the standard 
linear r-matrix algebra for a Hamiltonian system of natural form where the potential 



U is of second degree in the coordinates (see e.g. || |9], |T|, [L8|, |2T[). In recent years the 
study of completely integrable systems admitting a classical r-matrix Poisson structure 
with the r-matrix dependent on dynamical variables has attracted some attention |3], 
[8], jl3| . It is remarkable that the celebrated Calogero-Moser system, whose complete 
integrability was shown a number of years ago (c.f. [[HJ), was found only recently to 



possess a classical r-matrix of dynamical type ||. In this paper we study another 
example of a dynamical r-matrix structure. The potentials described are known as 
a parabolic family of potentials since they are permutationally symmetric potentials 
separable in generalized parabolic coordinates. They were introduced in |22|, and 
their connection with restricted coupled KdV flows was studied in [|T7| (see also [II for 



the isomorphism with KdV). However the r-matrix Poisson structure associated with 
higher degree potentials has not been discussed. The systems represent a generalization 
of a known hierarchy of two-particle systems with polynomial potentials separable in 



generalized (n- dimensional) parabolic coordinates (c.f. JTO, [TBI) 
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The polynomial second order spectral problem associated with the Lax representa- 
tion given below has been studied in a number of papers (c.f. f7|, |2~3}]). Here we reproduce 
all these results in the framework of a (2 x 2) Lax representation, which is a direct gen- 
eralization of that given by this being an effective technique to describe explicitly 
our class of integrable systems and to investigate their classical Poisson structure. We 
also note, that the Lax representation can be extracted from the results JTJ by some 
limiting procedure. We consider the system within the method of separation of vari- 
ables ITT] , |i~2| , |2ll which allows us to develop the classical theta-functional integration 
theory and to consider the associated quantum problem. The last problem is reduced 
to a set of multiparameter spectral problems which are a confluent form of ordinary 
differential equations of Fuchsian type. 

The central result of this letter is the description of the Poisson structure of the 
system by a dynamical linear r-matrix algebra. 



2 Lax Representation 



Consider the hierarchy of the Hamiltonian systems of n + 1 particles defined by the 
Hamiltonians 



H N (p 



n+1 



1; 



,p n+1 ; q 1 ..., q n+l ) = - ^ v\ + U N (q u . . . , q n +i) 
1 i=i 



(2.1) 



where the potentials fix the member of hierarchy by the recurrence relation 
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with the first trivial potentials given by Uq 



-2, = -2q n+1 - 2B, U 2 = -2q 2 n+1 - 



\ Y^l=\ <li- The general expression for the potentials Un at iV > 2 can be written in the 
form of a (N — 2) x (N — 2) determinant 



u 



N 



\N-l 



det 



/ In g-i go gi ■■■ gN-n \ 
/n-i —l g-i go ■ ■ ■ gN-e 
f N - 2 -l . . . g N - 7 



h 
V h 








-1 g-! 
... -1 / 



with 



g-i = Qn+i - B, g Ti 



Y,^Tql m = 0, 
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(2.3) 



f k = ^2Uig k -i-2, k = 3 
1=0 



,...,N, 



and Uo, U\, U2 as given above. The first nontrivial potentials are 
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2 



U, 
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(2.4) 
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+B Uq S n+1 - 2Bq 2 n+1 + q n+l £ gjj . (2.6) 



The potential (|2.5| ) is a many particle generalization of one of the known integrable 
cases of the Henon-Heiles system for which n = 1 (c.f. ||, |I4|). Analogously the potential 
( j2.6| ) is a many particle generalization of the system "(1 : 12 : 16)", known to be 
separable in parabolic coordinates (c.f. [|1(], |16[)- The system with the potential ( |2.5|) 



possesses two remarkable reductions: a) at q n+ \ = const, it reduces to the Neuman 
system which describes the motion of a particle over a sphere in the field of a second 
order potential and b) at q n+ \ = J2i=i qf it reduces to an anisotropic oscillator with a 
fourth order potential (see, for instance [0). The same reduction can be carried out 
for other members of the hierarchy. 

We look for a Lax representation in the form L N (z) = [M N (z), L N (z)], where 

L ~M = (S) -vt))' M " (z) = {qX) J)' (2 ' 7) 



with 



if ql 



U(z) = 4z-4q n+1 +4B-J2^TTi (2-8) 

i=i z + A 

V(z) = -lu(z) = 2p n+1 + ±^-, (2.9) 

W N (z) = -±U( z ) + U(z)Q N (z), (2.10) 

where Qn(z) is a polynomial of degree N — 2. The ansatz for the functions U(z), V(z), 
W(z) is a generalization of the corresponding ansatz constructed by Newell et al. to 
give the Lax representation for the integrable Henon-Heiles system. Here we introduce 
additional degrees of freedom n > 1, and higher degrees of the polynomial Qn are 
considered. See also |12| for the link to the su(l, 1)-Gaudin magnet which corresponds 
to a free n-dimensional particle separable in parabolic coordinates. 

It is possible to show that the Lax representation ( |2.7| ) is valid for all the hierarchy 
of Hamiltonian systems ( |2.1|) , ( |2.3j ) with the polynomial Qn( z ) and the function Wn(z) 



given by the recurrence relations 

Q N {z) = zQ N -i(z)-- ^ , (2.11) 

2 oq n+1 

W N (z) = W+(z)+W~(z), (2.12) 

W+(z) = zW+_ 1 -2U N ^ 1 , iV = 2,..., (2.13) 

W-(z) = E-^r, (2-14) 
where Un-i is the potential fixing the (N — l)-th member of the hierarchy and Q2 = 1. 
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We can easily solve ( 2.11 )-( 2.14j) to obtain formula for the functions Wn(z) and Qn(z) 



Qn(z) 
W N (z) 



JV-2 



N-l 



1 dU-N-k-l z k 

2 fc=0 dq n+ i 

JV-2 

- 2 £ u N -k-iz k + E 



pi 
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For example, for the first nontrivial cases we have 
Q3O) = z + 2g n+ i, 

n 

W 3 (z) = Az 2 + 4zq n+1 +4Bz + 4:q 2 n+1 + Y / Qt + Y, 



Pi 



k=i 



iZ + Ai 



for the many particle Henon-Heiles system and 



(2.15) 
(2.16) 

(2.17) 
(2.18) 



Qi{z) 



+ 2zq n+1 + 3ql +1 + - J2 Qi ~ 2Bq n+ll 



i=l 



W 4 (z) = Az 3 + ABz 2 +Az 2 q n+1 + Azq 2 n+1 + zJ2qi 



i=l 



+ 2q n+l Y,q--Y,Ml-ml + i + Y,-^T 

i=l i=l i=l Z + A i 



for the system with N = 4. 

The Lax representation yields the hyperelliptic curve C^ N ' = (w, z) 

Det (L^ N \z) — wl) = 



(2.19) 



(2.20) 



(2.21) 



generating the integrals of motion H^,F^\i = l,...,n. We have from ( |2.21|) and 
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16z N - 2 (z + B) 2 + 8H N + J2 



JV 
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N = 3,..., 



(2.22) 



where 
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(2.23) 



with lij = qiPj — qjPi, i,j = 1 n. The integrals of motion H N , F$ , 2 = 1 

are independent and have vanishing Poisson brackets. 



n 
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3 Separation of Variables 



To define the separation variables, i.e. the canonically conjugated variables 71*,^, i 
1, . . . , n + 1 and n+1 functions $j such that 



$ j (n i ,m,H N ,F8\...,FV>) = 



where H^, F$ are the integrals of motion in the involution, we use the scheme (c.f. 



12| , |19| ) P0| ) according to which these variables are defined in terms of the Lax matrix 
as 7Tj = V(fjLi), U(fii) = 0. The set of zeros = 1, . . .n + 1 of the function U(z) in 
the Lax representation ( |2.7| ) defines the parabolic coordinates given by the formulae 

n n+1 

q n+1 = ^2 Ai + B + ^2 fit, 

i=l t=l 



q m = -4 —fa / ; — , m=l,...,n. (3.1) 



The momentum 7T m canonically conjugated to fi m have the form 



n+1 



7r m = V(fii) = (i m n ^ m * m= 1, ...,n+ 1. (3.2) 

i=l,...,n+l 

The separation equations are of the form 

7Ti=w 2 (fii), £ = 1, ...,n + l, (3.3) 

where the function w 2 (z) is given by ( |2.22| ). 

The separation equations have two uses - to integrate the equations of motion in 
terms of theta functions and to quantize the systems. We mention here that canonical 
quantization in the space of separation variables leads to the following multiparameter 
spectral problem for the wave function of the system \1/ = TYj=i 

[-— + 16x N ~ 2 (x + Bf + 8A n+1 + £ A x . . . A n+1 ) = (3.4) 

2=1 £ 

with j = 1, . . . , n + 1 and the spectral parameters Ai, . . . , A n+ i The problem ( p.4|) has 
to be solved onn + 1 different intervals -"permitted zones". 



4 r— Matrix Representation 

It follows from the results of §2 and || that the system admits an r-matrix algebra. 
Let {-, -} be the standard Poisson bracket and {• ® •} be the standard Poisson bracket 
in the product of two linear spaces V 2 <S> V 2 . Then the classical Poisson structure for 
the hierarchy of dynamical systems described by the Lax operator L(z) with the entries 
(|2.8|) - (|2.10|) can be written in the form 

{L[ N) (x)®Lf\y)} = {r(x-y),L{ N \x) + L { 2 N \y)] 

+ [ S W(x,y),L{ N \x)-Li N \y)}, (4.1) 
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where L\ \x) = I ® L^ N \x), L\ ' = L^ N \x) ® I, J is the 2x2 unit matrix and the 
matrices r(x — y) and s^ N \x,y) are given by the formulae 



r(x — y) 



x-y 



-P, P 



^ 1 ^ 

10 

10 

V o o o i J 



y) 



2aN(x,y)S, S = er_ ® 





1 



(4.2) 
(4.3) 



with 



a N (x,y) 



Q N {x) -Q N {y) x 



N-l 



y 



N-l 



x-y 
1 N ' 3 



E 



y 



x-y 

N-k-l Q U 



N-k-1 



x-y 



dq, 



(4.4) 



n+l 



The equality (|4.1| ) contains all the information concerning the hierarchy of dynamical 
systems, in particular one can obtain a simple proof of the involutivity of the integrals 
of motion. The representation ( [2.7|) can also be derived from ( |4.1|) (c.f. |3|, |i~8|l). 

We write the relation ( f4.1| ) in the form 



{L?\x) <? Li N \y)} = [d^\x,y),L^\x)} - [4f (x,y), L™(y)] (4.5) 

Tij at i < j. 



with ^r^fU^^f 



The compatibility conditions for ( [4.5| ) (Yang-Baxter equations) have the form 

[4? (*, y) , 4? (*,*)] + [4? (*, y) , 4? (v, *)] + [4? y) , 4? (*,*)] 

+ {4*°(i/) ? - {Li N \z) ® 4? (*,!/)} 



+ [ C (a;,y,«) 1 Lj v) (y)-^(2 ! )]=0 



(4.6) 



and cyclic permutations. In this context L[ N \x) = L^- N \x) ® J 



L( N \y)®I,L 



I, Lf\y) 



I <S> I ® L^ N \z), and c(x, y, z) is some matrix dependent on the 



dynamical variables. 

If we denote Si 2 — S <S> I, £23 = I ® S, S13 = cr_ (g> / ® <j_ , where the matrix 5 is 
defined in ( (4.3| ). Then the following equality is valid for each member of the hierarchy 
of dynamical systems 

{L?\y)® S W(x,z)} - {Li N \ Z )®s£\x,y)} 

= 2/3 N {x,y,z)[P 23 ,S 13 + S12] 



dp N (x,y,z) 



[s,Lr\y)-Lf\z)} 



dq n +i 

with cyclic permutations. In (|4.7Q the matrix s = o_ <g> c_ <g> a_ and 

Q N (x)(y -z) + QN{y){z -x) + Q N (z)(x - y) 



f3 N (x,y,z) 



(x - y)(y - z)(z - x) 



(4.7) 
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Therefore the content of this letter can be interpreted as a finding a solution for the 
dynamical Yang-Baxter equation ( |4.6| ) which describes the evolution of the hierarchy of 
many particle one-dimensional Hamiltonian system separable in parabolic coordinates. 

Details of all these results will be published elsewhere ||. 

In conclusion we remark that the dependence of the r-matrix on dynamical variables 
can in principle be avoided by embedding the system into a completely integrable system 
with more degrees of freedom. This system can be described by the standard linear 
r-matrix algebra ( [4.1|) with s = 0. For example, for the case N = 3, set the functions 
U(z), V(z), W(z), to be the entries of the Lax operator L(z) 



U{z) = U(z), V(z) — V(z) — 8Vz — 8BV, 
W{z) = W(z) -lQV 2 x + Q-8Vp n+1 + 16V 2 q n+1 , 



(4.9) 



where the functions U(z), V(z), W(z) are given by the formulae (|2.8|) - (|2.10|) at N = 3 
and Q,V are new canonically conjugated coordinates with respect to the standard 
Poisson bracket. It is easy to see that the corresponding system satisfies the algebra 
( |4.1|) with the matrix s = 0. 

The corresponding enlarged dynamical system has integrals of motion Z 3 , 7i 3 , T\, % = 
1, . . . , n given by the formulae 



Is 

n 3 



(0 



m=l 



H 3 + 



X 3 (B-q n+1 )+V(Aj2Piqi + V^2q 2 



in n -i n in 

-is E £ + -p 2 (b E 9? - 1 E Ml - k+i E £) 

4 i=i i=i * i=i 8 i=i 



if n n ^ \ 

+ ( B - Qn+l) J^PiQi ~ E A iPi1i + jPn+l E ?i ' 

1 \ i=l i=l 4 i=l ) 



(4.10) 
(4.11) 



(4.12) 



where in ( |4.11| ),( ^.12| ) H 3 and are the integrals of motions of the many-particle 
Henon-Heiles system calculated by the formula ( |2.23| ). We can see that at V = 0, X 3 = 
the system reduces to the many-particle Henon-Heiles system. 

We expect that the analogous apparatus can also be developed for systems separable 
in elliptic coordinates, in which case we have to change the ansatz for U(z) to U(z) = 
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